New gamma ray signal from gravitationally boosted neutralinos at the galactic center 
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We discuss on the possibility that colliding dark matter particles in the form of neutralinos may 
be gravitationally boosted near the super-massive black hole at the galactic center so that they 
can have enough collision energy to annihilate into a stau pair. Since in some phenomenologically 
favored supersymmetric models the mass splitting between the neutralino and the lightest stau, one 
of the two scalar superpartners of the tau lepton, is a few GeVs, this channel may be allowed. In 
addition, staus can only decay into a tau lepton and another neutralino. We calculate the gamma-ray 
spectrum and flux generated by the tau pair discussing the observability of the obtained features. 
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Dark matter (DM) accounts for more than 80% of the 
mass of the Universe but its nature is still one of the open 
problems in Physics. In a widely accepted theoretical sce- 
nario, DM is formed by a weakly interacting massive par- 
ticle (WIMP) that has been in thermal equilibrium with 
Standard Model (SM) matter in the early Universe, leav- 
ing, after decoupling, the DM relic density as inferred by 
WMAP 0. In this light, supersymmetric (SUSY) exten- 
sions of the SM provide a natural WIMP candidate. In 
the minimal supersymmetric standard model (MSSM), 
R-parity conservation assures that if the lightest super- 
symmetric particle is the lightest of the four neutralino 
states-indicated as \ m t he following-, this particle is ab- 
solutely stable. In a phenomenologically favored scenario 
of the constrained MSSM (CMSSM), the stau coannihi- 
lation region (tcr) 0]j the lightest stau, fi, one of the 
scalar super-partners of the tau lepton, is close in mass 
to the neutralino. In the tcr parameter space the cross 
section for non-relativistic annihilation into fermions of 
the SM, XX ~ //) is typically small and results in a 
too large relic density. However, including the so-called 
coannihilation processes as for example %fi, T\T\ col- 
lisions, when the mass splittings of the involved particles 
are small, one can efficiently enhance the thermally aver- 
aged cross section (av), and, consequently, diminish the 
relic density to the measured value. 

The standard cosmological model predict that non- 
relativistic cold DM particles (v/c ~ 1CU 3 ) cluster into 
halos [i[ that contain baryonic matter. Since DM in the 
halo follows a certain mass distribution, the two-body an- 
nihilation processes can happen at a rate that is propor- 
tional to the DM mass density squared. Therefore, the 
highest chances to detect an observable indirect signal 
of their existence are attained in a region with high DM 
density, in particular, in the galactic center (GC). Among 



the various signatures from DM annihilation, gamma-ray 
signals have received much attention. A continuum spec- 
trum of secondary photons may arise from hadronization 
and decay of the annihilation products [5jand from ra- 
diation from final state charged particles Q. Direct an- 
nihilation into photons is also possible but only at loop 
level 0- 

The gravitational potential in the GC is dominated 
by a super-massive black hole (BH) with mass Mbh = 
4 x 10 6 M Q and Schwarzschild radius R s = 2GM BH /c 2 = 



4 x 10" 7 



pc 



Recently, the idea that a BH can act 
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as a particle accelerator has been proposed [9|. The 
highest center of mass frame (CMF) energies are ob- 
tained when the colliding particles approach the horizon 
on falling geodesies with opposite angular momentum per 
unit mass L/m x < L c = AG Mbh. /c, in the case of the 
Schwarzschild metric. The maximum possible value is 
yfa- 2V5m x [l,[lo| for a non-rotating BH, while it can 
be arbitrarily large for the Kerr BH In principle, due 
to this general relativity effect, new annihilation chan- 
nels into heavier states, kinematically forbidden for non 
relativistic particles, could be accessible. Additionally, 
a realistic calculation of an indirect DM signal in this 
scenario would also be determined by the particle escape 
function at distances close to the BH. For the horizon 
proximity this has been calculated in under the re- 
strictive assumptions of annihilation into two massless 
particles with isotropic angular distribution. 

In this work we show that if DM is formed by neu- 
tralinos with the characteristics of the tcr, a new domi- 
nant annihilation channel may be opened already for sub- 
relativistic neutralinos [l2[ boosted in the inner regions 
of the GC. 

Near the BH the DM density is described by a power- 
law p(r) oc r -7 , as we will discuss later. From the New- 
tonian approximation given by the Jeans equation, the 
root mean squared velocity is v(r) s» (GMbh/t) 1 / 2 f]~3l — 
[l5j . or in terms of the Schwarzschild radius, v(r)/c ~ 
(i?s/2r) 1 ^ 2 . Since a Keplerian orbit with L c would cross 
the horizon if the pericenter distance is less than r m j n = 
4i?s, hence we will consider safely r > 4i?s- In this way, 
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FIG. 1. Diagrams for stau pair production and decay in neu- 
tralino annihilation. 



for example, at r = 10i?s we have v/c ~ 1/V20 ~ 0.22. 
Requiring that the neutralino has the relic density in- 
ferred by WM AP, the relative mass splitting with the fi , 
S = (mf — m x )/m x , is typically less than 5%. In the 
CMF the energy threshold for stau pair production is 
s/s = 2E X > 2m?, that implies v/c > [1 - 1/(1 + <5) 2 ] 1/2 - 
With 5 — 2%, v/c> 0.197: there exists thus a range of 
radii where the kinetic energy is high enough to reach 
the threshold of the process XX t~t + . The maximum 



radius is given by r n 



1/2[1 - (m x /m f f]R s . 



If the mass splitting Am = m? 



is larger than 



the tau mass, m T — 1.777 GeV, the staus can only decay 
into the two body final state %r, see diagrams in Fig. Q] 
The neutralino is bino-like thus the vertices's x~X~Z and 
X~X~{h, H) are suppressed, while the vertex x _ ^i _r is 
not suppressed by mixing. In fact it is proportional to 
Z\\U\2, the product of the relevant neutralino and stau 
mixing matrix elements that are both close to one. The 
dominant diagrams in Fig. [T] are thus the ones with t, u 
channel exchange of the tau. At energies near the thresh- 
old the produced staus are slow thus the propagator 

-2E x E f + 2pf Px -m 2 T ) 
— m 2 ]. The cross sec- 



ru 



l/(p 2 x -p f ) 2 ~m 2 T = l/(m 2 f - 
is approximately l/[(rrif — 
tion, proportional to the square of this quantity is, thus, 
enhanced for mass splittings approaching the tau mass. 

We illustrate the above features in Fig. [5J where we 
show the relevant cross sections as a function of y/s, left 
panels, for four points of the tcr that are allowed by 
present phenomenological constraints. The values of the 
universal scalar mass m , gaugino mass m.1/2, trilinear 
scalar coupling Ao and the ratio of the two Higgs expec- 
tation values tan/3 that define the CMSSM parameter 
space are given in Table HJ The numerical computation 
was done using the interfaced codes MicrOMEGAs [l6l |. 
CalcHEP [l7[ and SOFTSUSY [3. The point A is similar 
to best-fit point found in [lj|, that predict a light Higgs 
around 119 GeV. Next we take three cases with Aq = 0, 
as usual for setting upper limits with LHC searches. In 
the point B the Higgs is slightly heavier. Points C and D 
have low tan/3 and the Higgs is around 115 GeV. In all 
cases, the cross section for stau pair-production clearly 
dominates by one or two orders of magnitude the cross 
sections for annihilation into fermions except when ^/s 
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FIG. 2. Annihilation cross sections in picobarn as a func- 
tion of the CMF energy (left panels) and annihilation cross 
section times the relative velocity in cm 3 /s as a function of 
the colliding neutralino velocity (right panels). The CMSSM 
points are specified in Table U and the annihilation channels 
are labeled in the upper left panel. 



corresponds to the heavy Higgs A, H masses where the 
s-channel propagators are resonant. In the right panels 
in Fig. [2] we show the annihilation cross section times 
the particle relative velocity, as a function the CMF ve- 
locity of colliding neutralinos; this is the quantity that 
enters in the calculation of indirect detection signals of 
the processes considered. Note that cr^rei for annihilation 
in staus, near the threshold, is at least an order of mag- 
nitude bigger than the freeze-out value 3 x 10~ 26 cm 3 /s 
and that these values corresponds to v/c ~ 0.1 — 0.2 that 
are just the ones that can be obtained with the gravita- 
tional boost discussed above. At low_tan /3, cases C and 
D, the "right" selectron and smuon (£r) tend to become 
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TABLE I. CMSSM points used in his work. The sign of /i is positive. The neutralino and stau masses are also reported. 





m (GeV) 


mx/2 (GeV) 


A (GeV) 


tan P 


m s (GeV) 


m fl (GeV) 


A 


452 


780 


1110 


41 


327.2 


333.6 


B 


858 


1780 





45 


789.0 


782.2 


C 


122 


600 





10 


247.6 


252.7 


D 


166 


805 





10 


337.3 


339.4 



degenerate in mass with ?i and the cross section for an- 
nihilation into pairs of these scalars is much larger than 
in the cases A and B. Although the masses of the parti- 
cles in point B are much heavier than in the other cases, 
the mass splitting is around 3 GeV and av is of the same 
magnitude. The same effect can be seen comparing case 
D with C. 

A possible signal of the opening of the new channel is 
given by the gamma-rays produced by the tau pair. The 
extension of the source is set by r max . This is too small to 
be resolved by present telescopes, thus we treat it as point 
source at the GC at a distance from us of D — 8 kpc. 
To evaluate the flux we first note that applying the small 
width approximation to the stau propagators, and given 
that BR(fi — ¥ r ± x) = 1, we have a(xX T ~ T+ XX) — 

&(xx Tif i)£i? 2 (fi -> tx) ^ cr(xx nn) = cff- 

We can thus evaluate the differential photon flux as 
^1 = M Tdrr*a-~(r)v dr)^-^(r) (1) 

In the integral we treat the distances in units of the 
Schwarzschild radius, thus r is dimensionless and a fac- 
tor i?| appears explicitly. We note some differences with 
the standard almost-static XX ~> t + t~ annihilation: (i) 
there is no factorization into a particle physics and as- 
trophysics factor because all the factors in the integrand 
depend on r through the velocity dependence. An in- 
tegration over the CMF scattering angle is implied in 
a?f(r) that is evaluated taking the exact spin averaged 
squared matrix elements from CalcHEP; (ii) we do not di- 
vide by 2 because the final state necessarily contains two 
neutralinos; (iii) the taus are not monochromatic and the 
spectrum changes with the collision energy y/s and ulti- 
mately with the distance, while in the static case the taus 
have an energy equal to the neutralino mass and the ra- 
diated photon spectrum is limited by _B™ ax = E T = m x . 

Before proceeding further we will discuss this last 
point. The taus energy spectrum can be easily obtained 
by applying a Lorentz transformation with parameters 
P = (1 — 4m?/s) 1 / 2 and 7 = ^/s/2mf, to the spectrum 
calculated in the rest frame of the stau. In this frame 
it has fixed energy and momentum, E* = (to? — to 2 + 
TO 2 .)/2TOf and p* = (E* 2 - to 2 .) 1 / 2 . The resulting en- 
ergy distribution is flat and limited, dN (fi — > x T )/dE T = 
1/AE, AE = £™ ax - E™ h \ with E™ in = y(E* - ftp*) < 
E T < £™ ax = j(E* + (3p*). The number of pho- 
tons with energy E 1 produced by a tau with energy E T 
is given by dN 1 /dx = 1/2 f(x) with x = E 1 /E T and 



f(x) = x~ 3/2 exp[g(x)} +glog[p(l - x)](x 2 - 2x + 2)/x. 
This formula was obtained in Ref. [2(} , to which we refer 
the reader for details, by fitting the photon yield from 
taus obtained with Monte Carlo simulations of the non 
relativistic process XX t + t~ . In this case the taus 
have equal energy, hence we use a factor 1/2 for the yield 
of one particle. The gamma spectrum at distance r is 
then obtained by integrating over the tau energy distri- 
bution, 

dN , . 1 f dE T ( E y \ „. , , . 

dE-^ = AW) J -t f {t) e{E ^ E ^ {2) 

E min (r) 

We have multiplied by 2 to obtain the yield of the pair. 
The Heaviside function takes into account that for fixed 
E 1 the integrand is zero if E 1 > E T . For this reason the 
photon energy cut-off is E™ n for each y/s. The absolute 
cut off when integrating over y/s will be at "/(r^^jE*. In 
fact asr-> r max , /3 -> and £™ in -> j(r max )E*. Note 
that 7(r max ) ~ 1 and E* ^ (mf — m x ), thus the cut off 
is indicative of the mass splitting between the neutralino 
and the stau. 

The last ingredient that we need to evaluate in Eq. ([1]) 
is the DM density profile. In [2l| it was shown that the 
adiabatic growth of the BH at the center of the halo 
causes a steepening, called spike, of the initial ha lo p ro- 
file toward the GC. Successive studies [H], [H, I22h28] | 
showed that considering physical effects such as scatter- 
ing of DM particles off stars, capture by the BH, self- 
annihilation and capture within stars during the evolu- 
tion of the DM distribution, results in a shallower profile, 
Psp(r) oc r~ 3 / 2 . For the SUSY models A-D, and gener- 
ally for all the tcr parameter space, the spin indepen- 
dent and spin dependent elastic neutralino-nucleon cross 
sections are in the range 10 -11 -10 -9 pb and 10 -8 -10 -7 
pb [29| , respectively, thus the energy lost by elastic colli- 
sions with baryonic matter is not likely to be important. 
At a certain distance from the GC the density reaches a 
value such that self-annihilation itself acts to stop further 
rising. It has been shown [27l HI] that this latter effect 
do not set the density to a constant value, usually called 
annihilation plateau or core, but results in a mild spike 
(MS) with Pms( r ) <x r~ x l 2 . In addition, the adiabatic 
compression of the gravitational potential caused by the 
baryons already in the bulge of the galaxy [3(| should 
also be taken into account. 

We hence model the profile considering that at the ra- 
dius r sp sa 0.2rh, with = 1.67 pc the influence radius 
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FIG. 3. Top panel: differential flux multiplied by Bot- 
tom panel: Differential flux (dots) and fit (dashed lines) with 
F{x) = F + Fi exp(- J B 7 / J Bi) + F 2 exp{-E y /E 2 ). F , F 1>2 , 
Ei : 2 are fit parameters. 



of the BH, DM density is given by a compressed Einasto 
profile, p sp — 5 x 10 6 GeV/cm 3 as in [22[ . From here 
the profile is given by p(r) = p sp (r/r sp )~ 7s P, 7 sp = 3/2, 
up to the radius r a where the density reaches the value 
Pa = m v /(av)otf. (cv)a is the annihilation cross section 
and tf = 10 Gyr 22, 23] is the elapsed time since the for- 
mation of the spike. Finally, the inner MS is p a (fAa)~ 7a , 
7 a = 1/2, up to the limit 4i?s- The radius r a is found 
by matching the two power-laws, r a = r sp (/3 a //j sp ) _1 / 7sp . 
The values of p a are between 10 11 — 10 12 GeV/cm 3 and 
r a are of the order ICT 4 pc for the CMSSM points A-D. 
In the cases C and D, av in the non relativistic limit is 
around 10 -29 cm 3 /s that would result in p a a factor 10 2 
lar ger than in A and B. Anyway, as was already noted 
in the cross section for annihilation into leptons is 
strongly velocity dependent as can be seen in Fig. [2] 



it rises rapidly reaching values around fO 27 cm 3 /s at 
v/c ~ 0.1 — 0.2, thus p a is of the same order as in A and 
B. 

In the top panel of Fig. [3] we show the differential flux 
multiplied by E^ to exhibit the behavior at the high- 
est energies near the cut off. In the bottom panel, the 
differential photon flux is given by the dots. We find 
that the spectral shape of the flux is well fitted with the 
sum of two exponentials, as shown by the dashed lines 
in the bottom panel. The functional form is F(x) — 
F a + Ft expi-Ey/Ex) + F 2 exp(~E 7 /E 2 ) with F , F h2 , 
-Ei, 2 fit parameters. 

The peculiar characteristics of the signal are: (i) its 
origin is in the innermost region around the BH where 
the DM distribution is given by the MS density; (ii) the 
differential flux presents a nearly exponential shape with 
a hard cut-off that is determined by the mass splitting be- 
tween the neutralino and the stau; (iii) the signal shows 
up at energies below 10 GeV. Interestingly, in this few 
GeV energy region there are some unexplained excesses 
over the known backgrounds [3l|, |32[ • However our pre- 
dicted signal is too feeble to account for them. It is ex- 
pected that after 5 years operation, the Fermi-LAT satel- 
lite reaches sensitivities of 10~ 10 — 0.5 x 10~ n photons 
cm~ 2 s _1 for energies between 0.5 GeV and 10 GeV [33| . 
The proposed signal can be one of the components ob- 
served by the collaboration. Furthermore, it might be 
discriminated by the new proposed experiment Gamma- 
Light [HI that should achieve a better energy and angu- 
lar resolution than Fermi-LAT in the interval (10 MeV-1 
GeV). 

In summary, we have shown that if DM is formed by 
neutralinos as described in the stau coannihilation re- 
gion of the CMSSM, stau pair production may be the 
dominant annihilation channel in the innermost region 
of the GC near the BH. We have further shown that 
the gamma-ray spectrum produced by the rr pair com- 
ing from the decay of the staus present peculiar features. 
This may further motivate the exploration of the GC by 
the Fermi-LAT satellite to achieve a better understand- 
ing of backgrounds and to look for a possible new indirect 
signal of the presence of a DM component. 
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